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Abstract. We present evidence for a strong correlation 
between the gaseous absorbing column density towards 
type 2 Seyfert nuclei and the presence of a stellar bar in 
their host galaxies. Strongly barred Seyfert 2 galaxies have 
an average N# that is two orders of magnitude higher than 
non-barred Sy2s. More than 80% of Compton thick (N# > 
~ 2 ) Seyfert 2s arc barred. This result indicates that 
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stellar bars are effective in driving gas in the vicinity of 
active nuclei. 

Key words: Galaxies: nuclei-Galaxies: Seyfert 



Large scale gravitational torques, such as bars and interac- 
tions, are thought to transport gas into the central region 
of galaxies and, specifically, in the vicinity of active galac- 
tic nuclei (Shlosman et al. 1990). Seyfert galaxies are the 
low luminosity subset of AGNs. According to the accreting 
supermassive black hole paradigm, the accretion rate in- 
ferred for Seyfert nuclei is low (10 _1 — lO _2 M0?/r _1 ) and, 
therefore, much fuelling from the host galaxy is not re- 
quired. Indeed, McLeod & Rieke (1995), Ho et al. (1997a), 
Mulchaey & Regan (1997) show that the occurrence of 
bars in Seyfert galaxies is not higher than in normal galax- 
ies. Yet, there is both theoretical and observational evi- 
dence that stellar bars do drive gas into the central re- 
gion of galaxies (Athanassoula 1992, Tacconi et al. 1997, 
Lainc et al. 1998). The resulting central concentration of 
gas might not be relevant for the fuelling process, but can 
play a role in the obscuration of the active nucleus, i.e. 
stellar bars could contribute to the obscuration that af- 
fects Seyfert 2s. This connection would be very important 
for the unified theories (Antonucci 1993). In this letter we 
tackle this issue by comparing the degree of obscuration 
in Sy2s with the strength of the stellar bar in their host 
galaxy. 



2. Dependence of the absorbing Njf on the bar 
strength 

Hard X-ray spectra can be regarded as the best tool to 
measure directly the absorbing column density in Seyfert 
galaxies. Recent surveys have significantly enlarged the 
sample of Sy2s for which an estimate of Njj is available, 
and have also reduced the bias against heavily obscured 
objects that plagued former studies (Maiolino et al. 1998a, 
Bassani et al. 1998, Risaliti et al. 1998). 

We restricted our study to obscured Seyferts in the 
Maiolino & Rieke (1995) sample, completed with 18 addi- 
tional objects discovered in the Ho et al. (1997b) survey. 
These two Seyfert samples are much less biased than oth- 
ers in terms of luminosity and obscuration of the active 
nucleus and, therefore, can be considered representative 
of the local population of Seyfert galaxies. 

For AGNs which are thin to Compton scattering (i.e. 
Nh < 10 24 cm -2 ) along our line of sight the N# can be 
derived from the photoelectric cutoff in the 2-10 keV spec- 
tral range, provided that the signal to noise is high enough. 
If the source is Compton thick then the direct component 
in the 2-10 keV range is completely suppressed and we 
can only observe the reflected component, generally lit- 
tle absorbed. As a consequence, Compton thick Sy2s are 
sometimes misidentified as "low absorption objects" when 
observed in the 2-10 keV range. However, the fact that the 



absorbing column density is actually higher (> 10 



21 



cm 



2 ) 



can be inferred from several indicators, such as the equiv- 
alent width of the Fe line and the spectral slope, and by 
comparing the X-ray flux with other isotropic indicators of 
the intrinsic luminosity. A more detailed discussion of this 
issue is given in Maiolino et al. (1998a) and Bassani et al. 
(1998). We just mention that sensitive spectra at higher 
energies, such as those obtained by BeppoSAX in the 10- 
200 keV band, can identify column densities in the range 



10 



25, 



1 (Matt et al. 1998, Cappi et al. 1998) or 
set a lower limit of 10 25 cra -2 (Maiolino et al. 1998a). 

In Tab. 1 we list all the Seyfert 2, 1.9 and 1.8 galaxies 
(i.e. those showing indication of obscuration) in the joint 
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Fig. 1. Distribution of the absorbing column density in 
Sy2s as a function of the bar strength 



ray data, and for which a bar classification is also available. 
With regard to the bar we generally favored the classifica- 
tion based on near-IR images (since these are less affected 
by extinction) and quantitative identifications of the stel- 
lar bars based on ellipticity, position angle, "boxy-shape" 
and brightness profile arguments as measured from digital 
(unsaturated) images. Otherwise the optical classification 
reported in the RC3 was adopted (de Vaucouleurs et al. 
1991). However, in almost all cases the RC3 identification 
of bars was in agreement with the other works. When- 
ever the RC3 classification was not in conflict with the 
near-IR and/or "quantitative" classification (i.e. in most 
cases), we further split barred systems in strongly barred 
(SB) and weakly barred (SAB) as reported in the RC3. 
Four objects not classified or classified as non-barred in 
the RC3, while their infrared images show indications of a 
bar; in these cases we adopt the classification SAB, with 
the exception of NGC1068, whose bar appears strong in 
the K-band image (Thatte et al. 1998). 

The completeness of this Sy2s sample is discussed in 
detail in Risaliti et al. (1998). We could not find obvious 



Table 1. Absorbing N# and bar classification for the Sy2s 
sample. 



Name Log(N H ) a Bar classification 5 
RC3 Othcrs c Adopted 



NGC1068 


>25 


SA 


bar 1 


SB 


NGC1365 


23.3 


SB 




SB 


NGC1386 


>24 


SB 


- 


SB 


NGC1808 


22.5 


SAB 


— 


SAB 


NGC2110 


22.5 


SAB 


- 


SAB 


NGC2273 


>25 


SB 


double-bar 2 


SB 


NGC2639 


23.6 


SA 


no-bar 3 


SA 


NGC2992 


21.8 


_ 


no-bar 5 


SA 


NGC3031 


20.9 


SA 


- 


SA 


NGC3079 


22.2 


SB 


- 


SB 


NGC3081 


23.8 


SAB 


double-bar 2 ' 7 


SAB 


NGC3147 


20.6 


SA 


- 


SA 


NGC3281 


23.9 


SAB 


- 


SAB 


NGC3393 


>25 


- 


bar 2 - 5 


SAB 


NGC4258 


23.2 


SAB 




SAB 


NGC4388 


23.6 


SA 


boxy-bulge 4 ' 5 


SAB 


NGC4507 


23.5 


SAB 


bar 2 


SAB 


NGC4565 


<21.8 


SA 




SA 


NGC4579 


20.6 


SAB 




SAB 


NGC4594 


21.7 


SA 




SA 


NGC4941 


23.6 


SB 




SB 


NGC4939 


>25 


SA 




SA 


NGC4945 


24.6 


SB 




SB 


NGC5005 


>24 


SAB 




SAB 


NGC5033 


20.9 


SA 


no-bar 3 


SA 


NGC5135 


>24 


SB 


bar 2 ' 3 


SB 


NGC5194 


23.7 


SA 




SA 


NGC5347 


>24 


SB 


bar 2 


SB 


NGC5506 


22.5 




„o-bar 3 


SA 


NGC5643 


23.3 


SAB 


bar 2 


SAB 


NGC5674 


22.8 


SAB 


bar 4 


SAB 


NGC7172 


22.9 




no-bar 3 


SA 


NGC7314 


22.1 


SAB 


bar 3 


SAB 


NGC7319 


23.5 


SAB 




SAB 


NGC7582 


23.1 


SB 




SB 


NGC7590 


<20.9 




„o-bar 3 


SA 


IC2560 


>24 


SB 




SB 


IC3639 


>25 


SB 


bar 2 


SB 


IC5135 


>24 




bar 2 ' 3 


SAB 


IC5063 


23.4 


SA 


„o-bar 2 


SA 


Circinus 


24.7 


SA 




SA 


IRAS0714 


>25 


SA 




SA 


IRAS1832-59 


22.1 


SA 




SA 


Mkl066 


>24 


SB 




SB 



a In units of cm -2 ; those values are from Bassani et al. (1998) and Risaliti et 
al. (1998). 

b SA = non-barred; SAB = weakly barred; SB = strongly barred. 

c References: 1-Thatte et al. (1997); 2-Mulchaey et al. (1997); 3-Hunt et al. 

(in prep.); 4-McLeod & Rickc (1995); 5-Alonso-Hcrrcro et al. (1998). 



there is no correlation between N# and the luminosity of 
active nuclei (Risaliti et al. 1998). 

Fig. 1 shows the distribution of N# for subsamples of 
increasing stellar bar strength, ordered from bottom to 
top. There is a clear tendency for N# to increase along 
the sequence. Tab. 2 quantifies this apparent trend. The 
median of the N# distribution^ increases by more than 
two orders of magnitude going from unbarred to strongly 
barred Sy2s. The confidence of the results is given by the 
Gehan test (Feigelson & Nelson 1985), that takes into ac- 
count also censored data: strongly barred and unbarred 
Sy2s have N# distributions that are different at a confi- 
dence level higher than 99%. Another impressive result is 
that more than 80% of Compton thick (N# > 10 24 cm~ 2 ) 
Sy2s are barred (13 out of 16), to be compared with ~ 55% 
of the general population and specifically of early type 
systems that typically host Seyfert activity (Sellwood & 
Wilkinson 1993, Ho et al. 1997a). Also 56% of Compton 
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Table 2. Properties of the N#-bar correlation in Sy2s 



Parameter Non-barred Weakly Strongly 

barred barred 



No. objects 


16 


16 


12 


Med. Log(Nj/) a 


22.1 


23.5 


24.4 


Compton thick: 








No. (%) 
Probability 6 


3 (18.7%) 


4 (25.0%) 


9 (56.2%) 




89% 


>99% 



a Median of Log(Nn) in units of cm~ 2 . 

b Probability for the N^r distribution of intermediate and 
strongly barred Sy2s to be different from the N_h distribution 
of non-barred Sy2s. 



thick Sy2s are hosted in strongly barred systems, to be 
compared with ~ 25% of the general population. 

3. Other indications 

Maiolino et al. (1997) found that Sy2s are characterized 
by a rate of non-axisymmetric potentials (including inter- 
actions and peculiar morphologies) about 20% higher than 
Syls, this difference appears significant. Hunt & Malkan 
(1998) find the occurence of bars in Sy2s not significantly 
higher than in Syls within the CfA and the 12/im sam- 
ples. In the samples of Ho et al. (1997a) and Mulchaey & 
Regan (1997) the occurrence of bars in type 2 Seyferts is 
10-20% higher than type 1 Seyferts. These results indicate 
that even if bars drive gas into the circumnuclear region, 
such gas does not reduce much the opening angle of the 
light cones. Yet, if the correlation between bar strength 
and amount of circumnuclear gas obtained for Sy2s ap- 
plies also to Syls, we would expect a large amount of gas 
in the circumnuclear region of barred Syls as well. 

The circumnuclear (cold) gas is expected to Compton- 
rcflcct the nuclear X-ray radiation. This Compton- 
rcflcctcd component should flatten the X-ray spectrum in 
the 10-30 keV spectral range. Therefore, within the bar- 
circumnuclear gas connection scenario depicted above, we 
would expect barred Syls to have a flatter spectrum in 
the 10-30 keV band. However, this test is subject to var- 
ious caveats. First, variability affects the slope of the ob- 
served spectrum because of the time-lag between the pri- 
mary and reprocessed radiation. Second, a fraction of the 
"cold" reflection is expected to come from the accretion 
disk. Third, the effect is expected to be small: the reflected 
component should contribute no more than ~30% in this 
spectral region. Fourth, to date spectra in this X-ray band 
are sparse and with low sensitivity. So far, the only (small) 
sample of Syls observed at energies higher than 10 keV is 
the one presented in Nandra & Pounds (1994), that use 
Ginga data. Their sample contains nine Syls whose host 
galaxy have a bar classification. As shown in Tab. 3 the 



hard X-ray spectra of barred Syls to be flatter than the 
unbarred Syls. 

Large amounts of circumnuclear gas in Syls could be 
detected via dust-reprocessed light in the infrared. More 
circumnuclear gas would imply more warm (AGN-heated) 
dust, hence more infrared emission relative to the intrinsic 
luminosity of the AGN (traced by the hard X-ray lumi- 
nosity). The mid-IR (~ 10/zm) is an excellent band to 
look for this excess, since the AGN IR emission peaks 
there (Maiolino et al. 1995). Also, by using narrow beam 
photometry it is possible to isolate the contribution of the 
AGN from the host galaxy. Within the scenario of the bar- 
circumnuclear gas connection, barred Syls are expected 
to show a mid-IR to X-ray flux ratio higher than non- 
barred Syls. Yet, several caveats affect this test as well. 
Both short and long term variability plague the reliability 
of the X-ray flux as a calibrator of the AGN luminosity. 
Equilibrium temperature arguments indicate that the dust 
emitting significantly at 10/xm should be located within 
the central 1-10 pc; therefore, excess of circumnuclear gas 
distributed over the 100 pc scale would not be probed by 
this indicator. Finally, even the small aperture (5") used 
in most of the groundbased mid-IR observations might in- 
clude the contribution from a central compact starburst. 
Tab. 3 reports the mean of the 10/xm/2-10keV luminos- 
ity ratio as a function of the bar strength for a sample of 
19 Syls. The 10/xm data are from Maiolino et al. (1995) 
and from Giuricin et al. (1995); the X-ray data are from 
Malaguti ct al. (1994), where we choose the datum closest 
in time to the mid-IR observation, to minimize long term 
variability effects. There is a tendency for barred systems 
to have a higher Lio^m/Lx ratio, though the spread is 
large and the statistics are poor. 

4. Discussion 

The important result of our study is that the absorbing 
column density of type 2 Seyferts strongly correlates with 
the presence of a stellar bar in their host galaxies. As dis- 
cussed in the Introduction, this result is not completely 
unexpected: stellar bars are very effective in driving gas 
into the central region, thus contributing to the obscu- 
ration of AGNs. On the other hand this gas should not 
play a major role in powering low luminosity AGNs, such 



Table 3. Properties of Syls as a function of the bar 
strength 



Parameter 


Non-barred 


Weakly 


Strongly 






barred 


barred 


(ri()_l8fe e v} 


1.69± 0.24 


1.31±0.10 


1.24±0.47 


(No. objects) 


(3) 


(3) 


(3) 


{log(L IR /L x )) a 


0.92± 0.07 


1.20±0.23 


1.58±0.80 



(No. objects) (5) (6) (8) 
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as Seyfert nuclei, given the lack of correlation between 
stellar bars and Seyfert activity (McLeod & Rieke 1995, 
Mulchaey & Regan 1997, Ho et al. 1997a). 

As discussed in the former section, the fraction of bars 
in Sy2s is only moderately higher (at most by 20%) than 
in Syls. So, the gas accumulated in the central region by 
the stellar bar increases the column density outside the 
light cones, but it docs not increase much the covering 
factor of the obscuring material. A possible explanation is 
that radiation and wind /jet pressure act to destroy or ex- 
pel molecular clouds that happen to enter the light cones, 
while outside the light cones self-shading allows molecular 
clouds to survive and pile up along our line of sight. An- 
other possibility is that the inflowing gas concentrates at 
dynamical resonances (eg. Lindblad resonances) forming 
obscuring tori. The solid angle that such tori subtend to 
the AGN depends on their inner radius and thickness, that 
in turn depend on the dynamical and kinematical proper- 
ties of the nuclear region, but are relatively independent 
of the amount of gas in the torus (if the gas self-gravity 
is significant the torus would flatten, thus actually reduc- 
ing its covering solid angle). On the other hand, the Nh 
through the torus depends linearly on its mass. 

Our result has also important implications on the 
scales over which the obscuring material is distributed. 
Large scale stellar bars can transport gas into the central 
few 100 pc, but they exert little influence on the gas dy- 
namics on smaller nuclear scales. As a consequence, a first 
interpretation of our bar-N/f link is that a large fraction 
of the obscuration in Sy2s occurs on the 100 pc scale. This 
is in line with results obtained from HST images, that as- 
cribe the obscuration of several Sy2s nuclei to dust lanes 
or disks a few 100 pc in size (eg. Malkan et al. 1998). Yet, 
one of the most interesting results of our study is that stel- 
lar bars appear very effective in making the obscuration of 
Sy2 nuclei so high to be Compton thick: more than 80% 
of Compton thick Sy2s are barred. As discussed in Risaliti 
et al. (1998), column densities larger than 10 24 cto~ 2 are 
unlikely to be distributed on the 100 pc scale, since the im- 
plied nuclear gas mass would exceed the total dynamical 
mass in the same region for several objects. This consid- 
eration generally constrains most of the Compton thick 
gas to be located within a few 10 pc from the nucleus. On 
the other hand, on these small scales the gas dynamics 
is expected to be little affected by the non-axisymmetric 
potential of a stellar bar in the host galaxy. The connec- 
tion between large stellar bars and Compton thickness of 
Sy2s thus requires some mechanism to link the dynam- 
ics on these different scales. Shlosman et al. (1989) pro- 
posed that if the large scale stellar bar collects in the cen- 
tral region a mass of gas that is a significant fraction of 
the dynamical mass, then the gaseous disk might become 
dynamically unstable and form a gaseous bar that could 
drive gas further into the nuclear region. A nuclear gaseous 



Seyfert 2 nucleus, in agreement with expectations. Ironi- 
cally, this galaxy has been classified as non-barred in the 
RC3 catalog; however the fact that it is edge on and lo- 
cated in the Galactic plane (where crowding and patchy 
extinction confuse the large scale morphology) might have 
prevented the identification of a large scale stellar bar. Al- 
ternatively, nested secondary stellar bars have also been 
observed in several galaxies and also in Seyfert galaxies 
(Mulchaey & Regan 1997), and are thought to be more 
stable (Friedli & Martinet 1993). An encouraging result 
in this direction is that the only two Sy2s showing ev- 
idence for an inner secondary stellar bar and for which 
the N# has been estimated are actually heavily obscured 
(see Tab. 1). However the opposite is not true: 5 other 
Compton thick Sy2s imaged in the near-IR by Mulchaey 
et al. (1997) do not show evidence for inner bars, though 
the limited angular resolution might have prevented their 
detection. 
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